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PREFACE 


The  work  described  in  this  publication  was  performed  by  the 
Mathematical  Analysis  Research  Corporation  (MARC)  under  contract  to 
the  Jet  Propulsion  Laboratory,  an  operating  division  of  the  California 
Institute  of  Technology.  This  activity  is  sponsored  by  the  Jet 
Propulsion  Laboratory  under  contract  NAS7-918,  RE182,  A187  with  the 
National  Aeronautics  and  Space  Administration,  for  the  United  States 
Army  Intelligence  Center  and  School. 

This  specific  work  was  performed  in  accordance  with  the  FY-87 
statement  of  work  (SOW  #2) . 
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2-D  vs  3“D:  Sine  of  Squared  Angular  Error  and  Weighted  Perpend icular 
Fixing 


INTRODUCTION 


J 


If 


1) 

2) 


two  dimensional  algorithms  work  on  variables  from  a  gnomonic  plane 
gnomonic  planes  are  based  on  the  same  point  of  tangency  both 
within  a  single  algorithm  and  between  the  two  and  three  dimensional 
algorithms 

then  it  is  reasonable  to  ask  the  question: 


How  much  difference  does  it  make  if  the  gnomonic  projections  are  taken 
before  or  after  fixing? 


This  question  will  be  put  into  mathematical  form  and  analyzed  in  this  memo. 
The  procedure  used  will  be  asymptotic  expansion,  on  the  assumption  that  error 
would  go  to  zero  as  the  radius  of  the  Earth  became  infinite  if  the  same 
gnomonic  projected  coordinates  were  Involved.  In  this  way  the  dependence  of 
the  difference  of  projections  on  the  scale  of  the  problem  may  be  seen.  This 
format  is  important  because  the  scale  restriction  is  a  significant  known 
parameter,  best  accounted  for  in  this  way. 

For  analysis  purposes  it  is  also  relevant  to  identify  how  error  models 
cross  to  the  gnomonic  plane.  They  also  may  or  may  not  correspond  well  with 
two  dimensional  error  models. 

The  two  fix  techniques  being  analyzed  are  the  sine  of  squared  angular 
error  and  its  approximation,  the  weighted  perpendicular.  The  weighted 
perpendicular  technique  is  used  by  the  Improved  Guardrail  V  (IGRV)  program, 
and  differs  from  the  sine  of  squared  angular  error  technique  in  interesting 
ways  that  will  be  discussed  later.  * 

These  topics  are  analyzed  in  detail  in  the  sections  that  follow. 


CONCLUSIONS 


Analysis  shows  there  is  little  difference  between  using  the  three 
dimensional  model  and  using  the  two  dimensional  model,  if  the  gnomonic  plane 
is  positioned  close  to  the  emitter.  By  little  difference  we  mean  that  the 
distance  between  the  fix  determined  by  these  two  methods  is  small. 

For  small  scale  problems  and  relatively  accurate  lines  of  bearing,  it  can 
be  shown  that  the  distance  between  the  fix  estimates  would  be  a  small 
percentage  of  the  maximum  miss  distance  of  an  LOB.  In  particular,  in  the 
appendix  it  is  shown  that  the  percentage  is  less  than  the  absolute  value  of 
the  maximum  of  the  miss  distance  to  the  radius  of  the  Earth,  ratio  squared, 
times  a  factor  to  account  for  the  size  of  the  base  relative  to  the  size  of  the 
problem.  Because  of  cancellation  owing  to  terms  with  different  signs,  the 
percentage  is  usually  less  than  the  squared  ratio. 
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The  small  difference  that  does  exist  between  the  fix  methodologies  only 
reflects  a  difference  in  the  weighting  of  the  different  data.  In  line  of 
sight  applications,  the  weighting  in  the  three  dimensional  method  makes  less 
sense  than  the  weighting  of  the  two  dimensional  method. 

METHOD  OF  ANALYSIS 

The  fix  methods  and  gnomonic  projection  are  not  affected  by  how  latitude 
and  longitude  coordinates  are  assigned  to  the  sphere  (to  within  spherical 
rotations).  Therefore  to  shorten  expressions,  we  will  start  by  assuming  that 
the  point  of  tangency  of  the  plane  with  which  the  gnomonic  projections  are 
taken  is  at  (0,0, R),  where  R  is  the  radius  of  the  Earth. 

z-R  is  the  gnomonic  Plane 

^Xi’^i*R)"the  sensor  location  in  the  gnomonic  plane,  (X^.Y^)  fixed. 
9i«the  angle  of  the  ith  line  of  bearing  in  the  gnomonic  plane 
(x,y,R)«the  location  estimate  in  the  gnomonic  plane. 

The  ith  plane  of  bearing  corresponding  to  0A  is  spanned  by  (X^.Y^.R)  and 

(sinei,cos9i ,0) ,  and  goes  through  (0,0,0). 

n  -  the  unit  normal  corresponding  to  the  ith  plane  of  bearing 
1  which  is  proportional  to  (X^ Y1,R)x(sin9i,cos9i,0) 

or  (-Rcos9i,Rsin9i ,Xicos0i-Yisin9i)//R2+(Xicc39i-Yisin9i) 2 

(x,y,R)(R//x2+yi*R2)-  point  on  the  Earth  with  gnomonic  projection  (x,y,R) 

(Xi,Yi,R)(R//X2+Y?>R2)-  point  on  the  Earth  with  gnomonic  projection  (X^Y^R) 

The  two  dimensional  likelihood  function  to  be  minimized  for  two  dimensional 
sine  of  angular  error  is 

L2  -  I[(x-Xi,y-Y1)(-cos9i,sin9i)]2/[(x-Xi)2+(y-Yi)2] 

The  three  dimensional  likelihood  function  to  be  minimized  (based  on  the 
inverse  image  of  the  two  dimensional  sensor  locations  and  bearings)  is 

L  -  E([(x,y,R)//x1+y2+R2  -( X  , Yi , R) //X* *Y2 *R2 ] • n> >2/Di 


Cl  -  ( -Rcos9 ^ ,  Rsin©^,  X^cose^-Y^ine^ 

d^  -  -(x-X^JcosQ^y-Y^sine^  -  signed  perpendicular  distance 
n>  -  normal  unit  vector  -  C^/||C^|| 

Fi  “  l|Ci||2  -  R2*(Xicos91-Y1sin91)2 

1/F1-  1/||C.|!2  -  1/R2-1/RI,[(Xicos91-Y1sin0i)2]  ♦  0(1/R6)  (as  R~>“) 


■eh 


Di  -  |  |[(x,y,R)//x2*y2  +  R2  -(X  .  ,  'i  i  ,  R)  //X2  +  Y|*R2  ]  |  |  2 
Ei  -  [(x,y,R)//x2*y2^R2-(Xi,Y. , R) //X2+Y2*R2 ]  C 


CALCULATIONS 


L  -  EES/CF.D^ 


n  min  L  -  E[-(x-Xi )cos9i+(y-Y . ) sin6i ]2/[ (x-Xi )2+(y-Yi )2]  -  min  L2 
—  >  00 


This  implies  that  the  limit  of  the  three  dimensional  fix  (as  the  relative  size 
of  the  Earth  to  the  scale  of  the  problem  goes  to  infinity)  is  the  two 
dimensional  equal  to  the  two  dimensional  fix. 

To  get  an  idea  of  the  actual  difference,  expand  the  Taylor  Series  for  the 
three  dimensional  fix  about  R«».  The  constant  term  is  the  two  dimensional 
fix.  The  first  derivative  term  shows  the  dependence  on  R. 

The  two  techniques  mentioned  above,  solve  the  problem  of  minimizing  L  and 
L2  differently.  The  sine  of  squared  angular  error  does  so  by  taking  the 
partial  derivatives  of  the  functions  and  solves  them  for  zero.  The  weighted 
perpendicular  technique  approachs  the  problem  by  taking  the  partial  derivative 
of  only  the  numerator  of  the  expression.  In  leaving  the  denominator  alone, 
this  technique  weights  the  value  of  the  derivative.  The  resulting  expression 
is  then  set  equal  to  zero  and  solved. 


3-D  FIX  METHOD 

We  will  need  to  use  the  correspondence  between  (x,y)  in  the  gnomonic 


projection  and  (x,y,R)//x2*y2+R2  In  the  three  dimensional  representation. 
If  x«f(l/R)  y«g(1/R)  denote  the  implicit  functions  defined  above,  then 


L  (f(1/R),g(1/R))  -  L  ( f ( 1 /R) , g( 1 /R) )  -  0 

x  y 


Hence  with  U-1  /R 


L  f’+L  g'+L  ,  -  L  f'*L  g  ’  ♦L  ■ 
xx  xy#  xU  xy  yyB  yU 


Thus, 


It  turns  out,  however,  that  LxU“LyU«0  at  R—  thus,  f '  (O)-g'  (0)-0. 


Therefore  to  get  the  first  order  term  we  need  to  compute  f"  and  g" .  Note 
first  that 


^’xx^"+‘Jxy8n+^xxx ^'^*2LXXyf'g'+LXyy(g')2*2Lxxuf'^2LXyug'*Lxu(j  -  0 


L  f  "  +  L  g"+L  (f’)2*2L  f'g’*L  ( g '  )2*2L  ,,f‘*2L  ,,g’*L  -  0 

xy  yy8  xxy  xyy  8  yyy  8  xyU  yyU8  yUU 
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However  since  f'(0)-g'(0),  this  reduces  to 


L  f "+L  g"+L 
xx  xy°  xUU 

L  f"-L  g" ♦L  „„ 
xy  yy  yUU 


0 
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L 
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f 

XX 

xy 

xUU 

L 

L 

r 

xy 

yy 

yUU 

at  the  point  of  interest  (U-0  or  R-«).  Thus  we  need  to  evaluate  this 
expression  and  take  the  limit  as  U->0  (R->»).  Since  we  are  examining  two 
techniques,  this  will  be  done  for  both.  However,  the  calculations  are  similar 
for  both  techniques,  so  the  majority  of  those  calculations  shown  will  be  for 
the  sine  of  squared  angular  error.  Wherever  a  difference  in  the  techniques  is 
of  interest,  the  calculations  for  both  will  be  shown. 


MATH  APPENDIX 


1  //R2  +a  -jEq  (-1  ^2)  r2(“1  /2-i)ai  .  (1/R)-(l/2)(1/R3)a+(3/8)0/R5)a2*  0(1/R?) 
Thus, 

(x,y,R)//x2>y2+R2-(x,y,R)[(1 /R)-(1 /R3) (x2+y2)/2+(3/8) (i /r5) (x2+y2)2+  0(1 /r7) ] 
(x,y,R)//xi+yi+R2-(Xi,yi,R)//X2+Y2+R2 

-(x-X.,y-Yit0)[1 /R]-([(x,y,R)(x2*y2)-(X1,Yi,R)(X2+Y2)]/2)[1  /R3] 

♦  (0(1 /R5) ,0(1 /R^) ,0(1 /R4 ) ) 

Thus  D.-[(x-Xi)2+(y-Yi)2][1 /r2] 

-C(x-X.  )(x(x2>y2)-X.(x2+Y2))*(y-Y.)(y(x2+y2)-Yi(X2+Y2))][l/R14] 

+  ([(x2  +  y2)-(x2  +  y2)]2/4)(;1/R4]  +  o(1/R6) 

And  so  1/D.-  [1/((x-X.)2+(y-Yi)2)][R2]  +  0(R4) 

Also, 

*L  -  [(x,y,R)//x2+y2+R2-(X1,Yi,R)//X2*Y2*R2]  C. 

-[ -(x-Xi)cos0i+(y-Yi) sine ^-[-xcosQ.+y sine i+X1cos9i-Yisinei ](x2+y2)/2[1 /R2] 

+  0(1/R4) 

-C-(x-X.  )cosei  +  (y-Yi)  sine  1]+[(x-Xi)  cose  ^(y-Yjsine^  (x2*y2)/2[l  /R2]-^  Od/R1*) 
-d  -d  (x2+y2)/2[1/R2]  ♦  0(1/R4) 

VlC<2E1,xDrEiDi,x)Ei/<FiDi):i  *  0(r2) 

VlU2Ei,yDi‘EiDi,y>Ei/(FiD!)3  *  0(r2) 

As  was  mentioned  earlier,  the  weighted  perpendicular  technique  (*)  utilizes 

a  different  L  and  L  .  Specifically,  L*  and  L*  do  not  have  D  and  D, 

x  y  x  y  i»x  i,y 

terms,  respectively.  Thus,  the  weighted  perpendicular  technique  has 

L‘-  e[(2E  E  )/(F  D  )]  ♦  0(R2) 
x  l,x  i  i  i 


L*-  E[  (2E  E  )  /( F  D  )  ]  ♦  0(R2) 
y  i » y  i  11 


where 


E  -  -cose  ♦[cose. ( x2*y2) *2x( (x-X. )cos9 ,-(y-Y  )sine  ) ]/2[ 1 /R2 ]  ♦  0(1/R4)  L 

1,  X  11  1111  * > 


-6- 


-  -cose . ♦[cosei(x2»y2)/2-xd . ][ 1 /R2]  ♦  OCl/R^) 

D  -  2(x-X. )[1 /R23"[x(x2+y2)-X  (X2*Y2)-(x-X  )(3x2*y2)+(y-Y  )2yx3[1/R^3 

1  ,  X  X  1X11  X 

+[(x2+y2)-(X2+Y2)]x[1/R^]  +  0(1/R6) 

-  2(x-x. )[1 /R2]-[ (x-X.) (X2+Y2)*(x-X. )(3x2+y2)+(y-Y. )2yx][1 /R^j-OC 1 /R6) 
E^  -  sine . +[-sine . (x2+y2) /2-yd^ ][ 1 /R2]  +  0(1 /R^) 

D.  -  2(y-Y.)[1/R2]-[y(x2  +  y2)-Yi(X2  +  Y2)  +  (y-Y.)(x2+3y2)  +  (x-X.)2yx][1/R1J3 
+[(x2+y2)-(X2+Y2)3y[l/R43  ♦  0(1/R6) 

2E  D.-E.D.  «[-2cose.((x-X  )2+(y-Y  )2)-2d  (x-X  )3  [1/R23 

ljXlllyX  X  1  X  X  X 

+2(cose.)[(x-X. )(x(x2+y2)-X. (X2+Y2))+(y-Y1)(y(x2+y2)-Y.(X2*Y2))3[1/Rl43 
-2 (cose.  )([(x2+y2)-(X2+Y2)32/10[1/Ri<3 
+[cose.(x2+y2)-2xd.3[(x-X. )2+(y-Yi)23[1/R43 
+d. (x2+y2)(x-X. )[1/R43 

+  d.[(x-X.  )(3x2*y2  +  X2>Y2)  +  (y-Y.  )2yx3[1/R1,3  ♦  0(1/R6) 

L  -Z[(2E.  D2E  -E2D.  D . +2E2  D2)~2(2E(  D.E  -E2D.  )D.  3/(F,d|)3 

XX  i,xx  111  i,xx  i  i,x  l  i,x  i  i  l  1,X  i,x  i  l 

L  -I[(2E.  D2E.-E2D.  D.+2E2  D?)-2(2E.  D.E.-E2D.  )D.  3/(F.d3)3 

yy  X.yy  1  1  i  i,yy  i  i.y  i  i,y  i  i  i  i.y  i.y  i  i 

L  -Z[(2E.  D2E.-E2D.  D.+2E.  D.  D.E.-2E.  D.  D.E +  2E.  E.  D2) 

xy  i,xy  ill  i,xy  l  i,x  i,y  l  l  l,y  i,x  i  x  i,x  i,y  l 

-2(2E.  D.E. -E?D .  )D.  3/(F.D?)3 

i,x  ill  i,x  i,y  i  l 

*'J'^2Ei,xyDiEi-EiDi,xyDi-2El  ,xDi  ,yDiEi-2Ei  ,yDi  ,xDiEi+2EiDi,xDi,y 

*2Ei,xEi.yD(]/(FlD?^ 

E.  «  E  -  E,  -  0  ♦  0(1/R2) 
i,xx  i , yy  i,xy 

D  -  2[1/R23  “  [(X2+Y2)  +  (3x2  +  y2)4,(x-X,)6x*(y-Y.)2y3[1  /R** 3  ♦  0(1/R6) 

X  f XX  XX  IX 

D.  *v“  -[(x-X.  )2y+2yx*(y-Y.)2x3[l/R1J3  ♦  0(1/R6) 
i » xy  l  i 

Da  -  2[1/R23  -  [(X2+Y2)*(x2^3y2)*(y-Yi)6y  +  (x-X1)2x3[1/R1<3  -  0(1/R6) 

Only  the  first  order  term  matters  as  R->»  for  these  terms.  Thus  at  R-*, 


E.  -  E .  -  E .  -0 

i,xx  i , yy  i,xy 


R2D ,  -  R2D  .  -2  R2D .  -0 

i,xx  i,yy  i,xy 


E,  --cose.  E.  -Sine,  R2D,  -2(x-X,)  R2D,  ,.-2(y-Y<) 


'iWw 


E1“[-(x-Xi )  cose i*(y-Yi) Sine  i]  R’D^C  (x-X1  )2+ (y-Y  .  )2] 


Ra/F  -1 


Lx,-ECdil-Di,,xDi*2D?.x,*‘,dloosei(t)lDl,x)'2COs29i  DI]/[RiDP 
-r[d2(6(x-X. )2-2(y-Y.)2)+8d. (x-X. )cosei((x-X. )2+(y-Yi)2) 

^cos2©^ (x-X.)2+(y-Y. )2)2]/[(x-X. )2*(y-Yi)2]3 

Lyy-r(d2(-DiyyD1+2D2jy)-i»disinei(DiDi>y)+2sin2eiD2)/(R*D3) 

-l[d2(6(y-Yi)2-2(x-X.)2)-8di(y-Yi)sine.((x-Xi)2^(y-Yi)2) 

*2sin2ei((x-X. )2+(y-Y. )2)2]/[(x-Xi)2+(y-Yi)2]3 

L  .  D.  D.  +  2d. (cose. D.  D.-sine.D.  D.)  -2sin0.cos0  D?)/(R*d3) 

xy  L '  1  1 ,  x  i  ,y  i  li,  y  l  ii,xi  i  ii  I 

-  rC8d2(x-X.  )(y-Y. )  ♦  Md^Cy-Y.  )cosei-(x-X.  )sinei)((x-X.)2  +  (y-Yi)2) 

-  2sin0iCose. ((x-Xi)2+(y-Y.)2)2]  /  ((x-X  )2*(y-Y  )2)3 


Thus, 

Exx  “  I[  (Mi  +  Ni  +  Oi)/(E2-*ri^)^] 

Lyy  '  E[<VVSl,/(cl‘nl)3:i 

Lxy  *  lC(W,i)/<t!',,|!,3] 

L  -  L 
yx  xy 

where 


Mi  -  d2(6e2-2n^) 

Ni  “  8diEicosei(E|+nf) 

0^  -  2cos20^(e2+n|)^ 

P1  -  d2(6r2-2e2) 

Q1  -  -8din1sin0i(e2+n2) 


-  2sin20^(e2+n|)2 

Ti  •  Sdfc  i*>i 


-  4d  .  (  n^ose^c^ine  . )  +  ) 

Vj  -  -2sine  .cose^e^  +  n2)2 

Using  similar  calculations,  the  weighted  perpendicular  technique  produces 

Lxx“  E(2^ef+h?)cos201  +  1Jcoseid1e  .  )/(e^  +  n^)2 

L*y-  Z(2(e2+n2)sin2ei  +  ifsineid .  n^/Ce^  +  n2)2 

L*  -  I(-2(e2  +  n2)sin6  .cose ,+^cose .d .n. )/(e2  +  n2)2 
xy  1  i  1  1  ill  i  i 

L*  -E(-2( e2+n? )sine  cose  -H sine .d.E)/(E?+n?)2 
yx  li  i  i  ill  li 


Returning  to  the  the  sine  of  squared  angular  error  technique,  recall  that 


L 

XX 

L 

yx 

-1  -i 

L 

yy 

-L 

xy 

L 

xy 

L 

yy 

Determinant 

-L 

yx 

L 

XX 

where  the  Determinant  -  L  L  -L  L 

xx  yy  xy  yx 


=ZI(M.Pj+M.QJ+M.SJ+N.Pj+N1Qj4N.Sj+0iPJ+0iQj+0iSj)/t(E2*n2)3(e2+n2)3] 


ZZ(TiTJ^2«TiUJ+2*T.Vj+2*UiVj+V.VJ+UiUJ)/[(E2  +  n2)3(E2  +  ri2)3] 


Terms  of  order  d^  may  be  ignored  for  relatively  accurate  LOBs.  All 
of  the  above  terms  are  at  least  order  d^  except  for  CLS^  and  V .  Thus, 

Determinant  -  [ ZEO^-V.Vj ) /[ ( e2+n2) 3( e2+n2) 333  4  0(di> 

-14ZZcos6isinej(e^*n^)2(e  j  +  rij  )2(coseisinej-sineicos6j  )/[  (e2  +  n^)3(E2+nj)3]+0(d 

-4zZcosei3inej3in(ej-ei)/[(E2  +  ri2)(E2  +  n2)  ]  +  0(di) 

-2ZZco36i3inej3in(0j-ei)/[ ( e2+n2 ) ( t2+n2) ] 

+2ZZcosej3in6  isin(  e^^-e  j  )/[(e2+n2)(e2  +  n2)]  4  0(di> 
-2ZZsin(ej-ei)(cose.sinej-cosejSinei)/[(E2*n2)(e2+n2)]  4  CKd^ 

F inally , 

Determlnant-2ZZsln2(ej-ei)/[(c2>n2)(E2+n2)]  ♦  0(d.) 


Lxx"[2Icos20i/(ci+ni)  +  0(di>3  4  0(1/R2) 


Lyy-[2Zsin261/(e2+n2)  4  0(d1)]  ♦  0(1/R2) 


Lxy-[-2Zsineicosei/(e2*nf)  4  0(dt)]  ♦  0(1/R2)  -  Lyx 

^  x-[-2coseidi/( e2+h2)  4  0(d2)]  ♦  0(1/R2) 


Li(y-[2sine1di/(E2+T12)  ♦  0{dJ)3  ♦  0(1/R2) 

With  these  assumptions( (d  ~0)and(R->») ) ,  the  IGRV  values  for  L*  ,  L*  , 

M  i  ^ ^  yy 


L*  L*  ,  L*  ,  and  the  Determinant*  coincide  with  those  above, 
xy’  i , x  i.y 


Computing  at  U-0  or  R-®  is  simplified  by  the  fact  that  the  terms 


making  up  L  can  be  chosen  so  that  L  «Ia  b  ...c  where  a  -b  «...-c  -0 

X  X  K  K  K  K  |  U  K  |  U  K  t  U 


at  U»0  (R-»).  Thus  mixed  terms  equal  zero  and  one  is  left  with 


LxUU“(Eakbk* ' ^k5 (ak,UU/ak  +  bk,UU/bk+*  * •+Ck,UU/ck) 


The  terms  used  in  the  two  kinds  of  products  making  up  L  ^  are: 


E.  ,  D .  /U2  ,  E.  ,  F  -U2  ,  and  U2/D.  . 

i  *  i,x  ’  i,x  i  i 


At  U*=0  or  R- 


Ei,uu/Ei-  -(x2^2) 


Ei ,  xUl/Ei  ,x 


-(x2+y2)+2xdi/cos6i 


(Di  x/U2)UU/(Di  x/U2)  “  -C(Xi+Yi>+3x2+y2]+2[(y-Yi)/(x-X.)]yx 


(Fi’U2)UU/(Fi’U2)  *  2(Xicosei-Yisinei)2 


(U2/Di)uu/(U2/Di)  - 

2C(x-Xi)(x(x2+y2)-Xi(X2+Y2))+(y-Yi)(y(x2+y2)-Yi(X2+Y2))]/[(x-xi)2+(y-Yi)2] 
-2([(x2^y2)-(X2  +  Y2)]2/l»)/[(x-Xi)2''(y-Y1)2] 

Ei.yuu/Ei,y  •  -<*2*y2)*2yV3ln6i 

(Di  y/U2)UU/(Di  y/U2)  “  "^Xf+Yf  )+3y2+x2]+2[  (x-X^/Cy-Y^  ]yx 


Terms  which  are  dependent  on  i  are  proportional  to  Lx  terms  and  hence  add  to 


zero.  Therefore,  at  U-0  (or  R«»)  and  d^_  0 


I(product)E.  ulJ/Ei  -  0 


Kproduct)  EliXUU/E1>x  -  0 


[(product)  E1>yUU/E1>y  '  0 


Thus, 


WEL1>IC<D1,X/U1W(D1.X/U')*(,/F-U’)UU/(1/F-U,)*(U'/D1W(U1/I>1)1 

LyUU-tLl,yl;(Dl.y/U')UU/<I>l.y/U',*<’/F'U>)UU/('/F-U,)*(U'/Dl>UU/lU‘/Dl)] 


In  order  to  calculate  f"  and  g",  where 
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L 

L 

-i  i 

XX 

xy 

xUU 

L 

L 

1 

xy 

yy 

yUU 

a  special  case  will  be  examined.  In  particular,  the  emitter  is  assumed  to  be 
at  (x,y)«(0,0).  This  is,  however,  a  logical  place  to  put  the  plane. 

Then , 

(Di,x/Ui)UU/(t,i,/U’)  -  -(Di,/U’V<Dl,y/Ul) 

('/Fi'u2)uu/(1/Fi*u2)  ’  "2di 

(UVDi)UU/(UJ/Di)  *  2[x2  +  y2-(x2+y2)/H]  -  6(x2*y2)M 
Therefore , 

LxUU-Zai,x('2d!  * 


Lyuu-I[Li.y("2<1!  *  (X!*,1)/2,J 


Since  the  weighted  perpendicular  technique  utilizes  a  different  L  and  L  , 

x  y 

it  should  produce  a  different  L,xUU  and  L  In  fact, 

Kw  1iX<-2<12  *  3/2(xf-Yf))] 

Lyw'  “Li,x(-2<1!  *  3/2(xfnf))] 

Notice  that  the  only  difference  in  these  values  and  the  ones  above,  is  that 
the  ( X2+y2 )  term  is  multiplied  by  3/2  rather  than  1/2.  This  is  due  to  the 
lack  of  a  (D^  )  term  in  the  equation  for  L*  (L*).  {Remember  that 

“’l.x^'V'Vx*’’’  - 

Finally, 

f"  -  [-1 /Determinant][L  L  MIt-L  L 

yy  xUU  xy  yUU 

g"  -  [-1  /Determinant][-L  L  ,.,+L  L 

yx  xUU  xx  yUU 

LyyLxUu‘LxyLyUU“C"II2sin2eicos6jdJ/(ei"ni)  ♦K2alne1cos019ineJdJ/(cf*nf ) 

+0(d2 ) ]  +  0(1/R2) 

-[2ZIsineisin(eJ-ei)dJ/(e2^n2)  ♦  0(d2)]  ♦  0(1/R2) 

~LyxLxUU+LxxLyUU’,t_II2sine1cos9icosejdj/(ef*nf)  ♦II2cos2019in0jdj/(ef*n?) 

+0( d2) ]  ♦  0(1/R2) 


-1 1- 


-[2ircoseisin(ej-ei)dj/(e^T,^)  ♦  0(d2)]  ♦  0(1/R2) 


The  only  difference  in  (-L  L  .  *L  L  )  and  (L  L  -L  L  )  and 

yx  xUU  xx  yUU  v  yy  xUU  xy  ynu'’ 

therefore  in  f"  and  g" ,  resulting  from  the  weighted  perpendicular  technique, 


is  that  both  are  increased  by  a  factor  of  three.  Thus, 


^yyLxUU-Lxy^yUu"^^sin®i5^n^0J-ei^dj'/^Ei4rii)  *  0<ai):l  4  0(,/r2> 
'LyxLxUu‘LxxLyUu‘tI£6coseiSIn(6j'6l,<1j/U?*',l)  *  0(di,:l  *  °(,'r2) 


Returning  to  the  sine  of  squared  angular  error  technique,  it  follows  that  f" 
and  g"  can  be  expressed  as 

f"-  -[2ZZsin0.sin(ej-e.)dJ/(e2+n2)]/{2lj:sin2(ej-ei)/[(e2+n2)(e2+n2)])  ♦  0(di) 

-  -[ZZsineisin(ej-61)dJ/(e2+n2)3/{ZZsin2(ej-ei)/C(e2*n2)(e2+n2)]}  *  CKd^ 
g"-  -[2ZZcosei3in(e  -6  )d  /(c2*n2)]/{2ZZsin2(e  -e.)/[(e2+n2)(e?+n2)]}  ♦  0(d) 

-  -[ZZcose i3in( e^-ei )d^/(e2+n2) 1/ ( IZsin2(0j-0j)/[ ( e2+n2) (e2*n2) ] }  ♦  0(d1> 


For  the  weighted  perpendicular  technique  used  in  the  IGRV  code,  a  different  f" 
and  g"  are  produced. 

f"»  ~3l  ZZsine^sinC 9j-ei )d^/( e2-*-n2)  ]/(  ZZsin2( e^-e^ ) /[  ( e2+ n2)  ( e2  +  n2)  ] )  ♦  (Kd^) 
g"-  _3[EIcoseisin(ej-ei)dj/(e2+n^)]/{ZZsin2(6j-ei)/[(e2+n2)(e2+nj)]}  ♦  CKd^ 

Since  the  only  difference  in  the  terms  for  f"  and  g"  between  the  two 
techniques  is  a  constant  multiplier,  in  simplifying  the  two  terms,  we  need 
only  concentrate  on  the  sine  of  squared  angular  error  technique. 

Note  that  if  i  0  for  all  i,  then 

minfa^/b^  i  Za^/Zb^  l  maxfa^b^} 
i  i 

Additionally, 

-sin(ej-ei)  -  sin(ei-6J) 

| | (sine^.cose^) | |  -  sin201  *  cos20i  -  1 


In  order  to  investigate  the  size  of  the  f"  and  g"  terms,  we  study  the  norm  of 
(f",g").  Using  the  facts  mentioned  above,  we  can  do  the  following. 


||(f",g")||  S  1z!j1|  |  (sln©1,cosei)  |  { 


j  £?[  sin(e  ©  j  )  d  j  /  ( c  **nj)  ] 
jlij’(sin2(el-ej)/[(eJ  +  np(ej*nj)]) 


£  max  jE^sin^-ejJdj  /  ^"[sin^e^e^/U  H,5| 


Now,  fix  the  i  such  that  the  expression  is  maximized.  Thus, 

I  |(f".g")|  I  *  |ji!^sin(ei-ej)dJ  /  jz"[sin2(ei-eJ)/(E2+n2)]| 

£  max{dj}jl!j1sin(ei-ej)  /  { [  1  /max{c2+nj  }  ]j  E^sin2  ( 6^6^ ) } 

j  J 

£  |  C  J  r*^sin(  0  e  J  )  /jE!j1sin2(ei-ej)3  max{dj}  max{e2+n2)j 

J  J 

Furthermore  since  nEX2  i  C ZXi l2 ,  calculations  show  that 
IX./ZX2  S  n/EXi  -  1/X 

Therefore, 

||(f".g")||  £  |max{d^.}max{e2+n2}  /  avgCsinCej-Bj ) ] 

J  J 


Finally,  the  distance  between  the  fix  produced  by  the  two  methodologies  (two 
dimensional  versus  three  dimensional)  is 

distance  S  | |(f",g")| |2  [1/R2] 

Si  max{dj}max{Ej  +  nj  }  /  avg[sin(6i-6j )  ]  1 2[1 /R2] 


w 


